Flexuron, a self-trapped state of electron in crystalline membranes 
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Self-trapping of an electron due to its interaction with bending fluctuations in a flexible crystalline 
membrane is considered. Due to the dependence of the electron energy on the corrugations of the 
membrane, the electron can create around itself an anomalously flat (or anomalously corrugated, 
depending on the sign of the interaction constant) region and be confined there. Using the Feynman 
path integral approach, the autolocalization energy and the size of the self-trapped state (flexuron) 
are estimated. It is shown that typically the size of the flexuron is of the order of the wavelength 
of fluctuations at the border between harmonic and anharmonic regimes. The flexuron states are 
connected with the fluctuation tail of the electron density-of-states, the asymptotic behavior of this 
tail being determined by the exponent of the renormalized bending rigidity. 

PACS numbers: 71.23.An, 73.22.Pr, 63.20.Ry, 68.60.Dv 
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I. INTRODUCTION 

Studies of statistical mechanics in two dimensions 1 
have been strongly stimulated recently by the discov- 
ery of graphene^, the first truly two-dimensional crys- 
tal and the simplest possible membrane. Experimental 3 
and theoretical 4 - demonstrations of intrinsic ripples on 
graphene due to thermal bending fluctuations have ini- 
tiated numerous works on the effect of the ripples on 
the electronic properties of graphene^— due to coupling 
between electronic and lattice degrees of freedom. This 
coupling originates from several factors, such as the de- 
pendence of electron hopping parameters on interatomic 
distances and on the angles between chemical bonds as 
well as from the redistribution of electron density in a 
deformed membrane (see for review Refs. iEil). 

However, graphene is just the first representative 
of a huge class of two-dimensional crystals includ- 
ing broad-gap semiconductors such as hexagonal boron 
nitride^, graphane (hydrogenated graphene)« and 
fluorographen e 14 ' 15 . Electrons in such materials are just 
conventional nonrelativistic quantum particles and not 
chiral Dirac fermions like in graphene, and one cannot 
expect there manifestations of exotic effects of corruga- 
tions such as gauge fields 9 . At the same time, another 
interesting physics specific of the two-dimensional case 
arises which will be the subject of the present work. 

It is known since many years that the interaction of 
a charge carrier in a semiconductor with some order- 
parameter fluctuations can drastically change its state 
leading to self-trapping, or autolocalization^—. This 
phenomenon is of crucial importance, for example, for the 
phase separation in magnetic semiconductors and colos- 
sal magnetoresistance materials 2 ^— , where the magneti- 
zation plays the role of order parameter. Since the band 
motion of the electron is easier (and, hence, the band- 
width is larger) for a ferromagnetically ordered state the 
electron in antiferromagnetic or magnetically disordered 
surroundings creates a ferromagnetic region and turns 
out to be self-trapped in this region (different names are 



used for this state, such as spin polaroid, ferronii, or 
fluctuoni 7 -). If concentration of electrons is large enough 
it leads to a formation of ferromagnetic state in the 
whole crystal, via so called double exchange, otherwise 
the phase separation happens, all electrons being trapped 
in the ferromagnetic regions. The crucial mechanism is 
the dependence of electron hopping parameter on the an- 
gle between magnetic moment on neighboring site o 16 ' 20 . 
Importantly, the critical point where fluctuations of the 
order parameter are the strongest is the most favorable 
for these "fluctuon" effects^^^^. 



Since the hopping parameter in a fluctuating mem- 
brane depends on the angle between normals^ one can 
expect a similar physics. If the electron motion is the eas- 
iest in the flat membrane, self trapping in an anomalously 
flat region is possible, in the opposite case an anoma- 
lously crumpled region can arise. Since membranes are 
characterized by very strong bending fluctuations and 
very slow, power-law decay of the corresponding cor- 
relation functions^ the situation is formally similar to 
that for autolocalization in the critical point in three- 
dimensional systems. I will study here this autolocal- 
ized state which can be called "flexuron", by analogy 
with "ferron", or "fluctuon". Technically, I will use the 
path integral approach 2 ^— as was applied to the fluctuon 
problem in our previous work s 18 ' 19 ' 23 ' 24 . This method is 
an adaptation of the seminal work by Feynman on the 
theory of polaroid. 



It will be shown that the flexuron forms in crystalline 
membranes at a large enough coupling constant and its 
size is of the order of the length scale where interactions 
between bending and stretching modes become impor- 
tant. The number of electron states in the fluctuation 
tail of the electron density of states where all states are 
(auto)localized, even without extrinsic disorder, turns 
out to be proportional to the temperature. 



II. DESCRIPTION OF THE MODEL AND 
MEAN-FIELD ANALYSIS 

In the continuum medium theory^ the membrane is 
described by the out-of-plane-deformation field h(x,y). 
The unit normal vector at a given point is 



n(x,y) 



fdh _dh ^ 



i + \vhY 



(i) 



where V is the two-dimensional gradient. In the ground 
state, the membrane is supposed to be flat, and all nor- 
mals are along z direction (Fig. la). At finite tempera- 
tures, the membrane is fluctuating, and the normals are 
no more parallel (Fig. lb). Let us assume, for simplic- 
ity, that the electron spectrum is isotropic near the band 
minimum. Then, by symmetry, the shift of the band 
edge due to bending fluctuations should be proportional 
to |V/i| 2 and the electron Hamiltonian takes the form 



ip(r) 



2m 
\Vh(r) 



- 9V M , 



(2) 



where r = (x, y), m is the effective mass, and g is the cou- 
pling constant, of the order of the electron bandwidth. It 
can be both positive or negative, as a result of the inter- 
play of different contributions. If the band is narrowed 
under crumpling, g > 0; the edge of the spectrum in this 
case corresponds to a homogeneous state with <p = 0, the 
energy of this state is E = 0. We will start from the case 
of positive g; for negative g the Hamiltonian ([2]) should be 
modified, to bound the ground-state energy from below. 
This case will be considered later, in Section 4. 

In the mean-field approximation, at finite tempera- 
tures the band edge is shifted up by 



E = g(\VK 



(3) 



The inclusion of fluctuations, however, changes this con- 
clusion dramatically. For the case of spin polarons, it was 
demonstrated already by Brinkman and Ricei£ that the 
band edges do not depend on the degree of spin disorder. 
The band is narrowed in antiferromagnetic or magneti- 
cally disordered state in comparison with the ferromag- 
netic one but not due to the shift of the band edges, just 
the density of states (DOS) becomes larger at the mid- 
dle of the band and exponentially small near the edges, 
similar to the so called Lifshitz tails 29 in disordered sys- 
tems. One can expect the same behavior for our case 
(Fig. lc). The mean-field energy Q is, rather, an inflec- 
tion point of DOS; the states between E = Eq and E = 
are associated with rare events, namely, correlated fluc- 
tuations such that in some large areas the membrane is 
flat and normals are parallel. The tail of the DOS below 
Eq is associated with flexurons. The autolocalization re- 
gion becomes larger when the energy goes closer to the 
band edge E = (and the edge itself corresponds to the 
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DOS(E) 
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FIG. 1: Ground state of the membrane, all normals (black 
arrows) are along the z axis (a), and a snapshot of a fluctuat- 
ing state at finite temperatures (b). Panel (c) shows a sketch 
of the density of states (DOS) at zero temperature (dashed 
line), at finite temperatures within mean-field approximation 
(solid line) and the fluctuation- induced tail (dotted line). We 
assume g > 0. 



completely flat state of the whole membrane, cf. Refs. 

mil). 

In crystalline membranes, the anharmonic effects are 
essential, namely, the coupling between bending and 
stretching modes; without this coupling the membrane 
turns out to be crumpled at any finite temperatures^. In 
the harmonic approximation which is valid for not too 
small wavevectors q of fluctuations, the Fourier compo- 
nent of the normal-normal correlation function G (q) = 



3 



\hq\ ) at temperature T is given by the expression^ 



G(q) 



T 
Kq 2 



(4) 



where k is the bending rigidity. The use of this expres- 
sion to calculate ^|V/i| 2 ^) leads to a logarithmically di- 
vergent result. The point is that the expression (J4j) is not 
applicable for small enough q where the anharmonic ef- 
fects become dominant. The coupling between the bend- 
ing and the stretching phonons in crystalline membranes 
leads to the renormalization of the bending rigidity and 
its growth with q decrease. Perturbation analysis shows 
that the anharmonic effects become dominant at 



q<q* = V/3TF/87TK 2 (5) 
where Y is the Young modulus ("Ginzburg criterion"—^). 



3 



For smaller q, 



G(q) 



.4 



where qo = \fYj k is of the order of inverse inter- 
atomic distance and 77 is the exponent of the renor- 
malized bending rigidity. Within the self-consistent 
screening approximation 3 -^— r\ ~ 0.821. Recent Monte 
Carlo simulations 33 and functional renormalization group 
analysi a 34 ! 35 give r\ ~ 0.85. 

The factor A can be found from the matching of ex- 
pressions (HJ and ([6]) at q = q*: 



where A = In {qo/q*) ~ (l/2)ln(«/T) originates from 
the cutoff at small q and we assume T< k (for example, 
^ for graphene k ~ let/ 4 ). Thus, within the mean-field 
approximation, the shift of the band edge is proportional 
to temperature. Physically this effect is observable (i.e. 
not hidden by thermal smearing) if |2?o| > T, that is, 
\g\ A k. Further we will assume that this condition is 
fulfilled, e.g., due to a large logarithm. 



A 



T 



1-17/2 



(7) 



where a is a numerical factor. 

With logarithmic accuracy, one finds 



Eq 



TgA 
2ttk 



III. FLUCTUATION EFFECTS IN THE CASE 
OF ELECTRON PREFERRING FLAT 
SURROUNDING 



To take into account the effect of fluctuations on the 
electron energy spectrum I will use the path integral 
formalisms^, following RefeIHIll2Hl 

The partition function of the whole system (the elec- 
(8) tron plus the fluctuating field tp) is represented as 



Z = Tre- pH f(!P)-PH.<r,<p) _ Zf ( Tr r T T exp 



H e (r,<p (t,t)) <h 



(9) 



r 



where j3 = T 1 is the inverse temperature, H e is the 
Hamiltonian ©, Zf = Tr v e" /3W ^ ¥ ' ) is the partition func- 
tion of the field, Ht is the corresponding Hamiltonian, 
tp ( r ,r) = e TU fM<p(r) e- TH s(v) anc j 

(A(v)) f = ±Tr <p e-^fMA(<p) (10) 
z f 

is the average over the field states. Using the Feynman 
path-integral approach 25-28 and taking average over tp 
yields for the electron-only free energy 



P P J r (0)=r(P) 

where the effective action of electron is 

5[r(r)] =5 [r(r)]+5 mt [r (r)] 

and 

P 



; 2?[r(r) 



exp (-S int [r (t)]) = ( exp 



P 

g / p[r(r)]dT 



(11) 
(12) 

(13) 
(14) 



(we put temporarily h = m = 1). The average can be 
formally written as a cumulant expansion 



Smt — — y, — r 

* — ' ml 



P 



lC m (r (n) ...r (r m )) dri...dr m 
(15) 

where /C TO (ri; ...;r m ) are the m-th cumulant correlators, 
defined recursively by 

Jd (r 1 ) = (^(r 1 )) / , 
/C 2 (r i; r 2 ) = (n) ^ (r 2 )) / - /d (r x ) K x (r a ) , ... (16) 

etc. 

To estimate the electron free energy J 7 we use the same 
trial action as in Refs. 18. 19.23, <S t — So + S pot where 



Spot - 4/? 



rP rP 

/ / [v{r)-v{cj)] 2 dTda, 
Jo Jo 



(17) 



the oscillator frequency ui being trial parameter (in con- 
trast with 

RefSlH we do not introduce any retardation 
in the trial action since it is important only for the case 
of quantum fluctuations of the field tp which is beyond 
the scope of the present work). It describes in a trans- 
lationally invariant way a self-trapped state if /3huj 1, 
in this case the size of the self-trapping region is of the 
order of the zero-point oscillations of a harmonic oscil- 
lator, I = \J h/2muj. The trial action St gives an upper 
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estimate of free energy for our problem. The Peierls- 
Feynman-Bogoliubov inequality reads 



1 



^int Spot) i 



(18) 



where J- t is the free energy corresponding to the trial 
action St, which is equivalent to 



F<Ft-— / ([r(r)-r(a)] 2 ) drda 
Jo Jo x 1 4 



U! 

4/3 2 



m=l J JO , =1 



(19) 



To proceed, we will pass to the Fourier transforms of 
the cumulants /C m (qi, .., q m _i) and take into account 
that for the Gaussian trial action St one has 



exp 



•E 



qjq/c 



([r^O-r^)] [r (r fe ) - r (r m )]) t 



(20) 



For the autolocalized states the variational parameter 
w satisfies the inequalities (3lj 3> 1 and ui <C where 
V4 7 is of the order of the electron bandwidth (the last 
inequality just means that the flexuron size I = \j\f2uj 
is much larger than the interatomic distance, otherwise 
our continuum-medium description is not applicable). 

Due to the first inequality one hasi£ 



(r { n )v{T k )) t ^2l 2 5 lk 



(21) 



and Eq.(|20p is simplified: 

T -\~ £ ^(.9/3) m E^m(qi,...,q m -i)x 

m— 1 



exp 




(22) 



In harmonic approximation, the second cumulant can 
be calculated using Wick's theorem: 



JC 2 (q)=^G(q)G(p-q) 



(23) 



Within the self-consistent screening approximation^ (ne- 
glecting vertex corrections) this expression is also sup- 
posed to be correct. Note that this approximation is in 
reasonable agreement with the results of Monte Carlo 



simulations^. We will use it here. Substituting Eqs.Q 
and © into Ea.(|23"|) one finds 



!C 2 (q) 



1 T\ 2 \n{q/q ) + 

tt UJ 9 2 , q>q 



q^q 2 - 



where a' is a numerical factor. 

Let us start estimating Eq. (l22"]) in the Gaussian ap- 
proximation, that is, taking into account only terms with 
772 = 1,2, then 

^<£o + -^-^$> 2 (q)expH Z 2 / 2 ) (25) 



Now we have to minimize the right-hand-side of Eq. (|25[) 
to find the size of flexuron Iq and the autolocalization 
energy. 

Let us assume first that q*l <C 1. In this case one can 
rewrite Eg. (|2"5]l as 



3TYH 2 



Tg 2 . 2 

■ In x, 



32tt 2 k 2 



(26) 



where x = (q*l) and we have restored the missing con- 
stants h and 771. 

The last two terms in the right- hand- side of Eq. (|2"6")) 
are proportional to temperature. The minimum of this 
function exists if 



g 2 m e 
V= 3^Y¥ > 2 



(27) 



(e ~ 2.71828...) which is a criterion of autolocalization 
in this approximation. Of course it is valid only with 
an accuracy of some numerical factor since we kept only 
leading logarithms in our estimations. It is hardly to 
expect that v is much larger than one, thus, an optimal 
value of x is of the order of unity, beyond the formal 
limit of the approximation (|2T)1) . The analysis from the 
opposite limit, q*l ^> 1 gives the same result which is 
not surprising keeping in mind that the expressions (p4|) 
match at q w q* . Therefore, 



lo * I/*?* oc T-V2, 
Topt — Eq — const ■ T 



(28) 



This is the main result of our consideration, and it is 
amazingly simple: the size of the self-trapping region is 
of the order of 1/q* (which is much larger than the inter- 
atomic distance, so the continuum medium description 
works). Thus, inhomogeneities induced by electrons in 
fluctuating membrane should have the same spatial scale 
as the crossover length from the harmonic to the anhar- 
monic regime. 

The higher-order terms in the cumulant expansion in 
Eq. p2|) can be estimated from the expressions in the har- 
monic regime. Each next couple of the Green functions 
G (q) gives an additional factor 1/q 2 , together with ad- 
ditional integration over q (which can affect the powers 
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of logarithms) and, importantly, the additional small pa- 
rameter (T/k) , which is only partially compensated by 
the additional factor fig. Thus, the cumulant expansion 
has a small parameter Tg/n 2 and the results of the Gaus- 
sian approximation are reliable. 

A compact analytical expression for the density of 
states N(E) can be obtained for the Lifshitz tail regime, 
that is, very close to the true edge of the spectrum E = 0. 
To this aim, one can use the method of Refill For the 
range of energies E 3> T ', relevant flat regions have a 
size much larger than 1/q*, thus, we are in the strongly 
anharmonic region. From scaling considerations^, one 
can postulate that for all m, as well as for m = 2 (cf. 
Eq.flU), 

K m (qi, q m -i) = a (2_2,?)(m_:L) /C m (aqi, aq m _i) . 

(29) 

Further consideration just repeats that of Refl23l for the 
case of XY-model. The result for the Lifshitz tail of the 
density of states is 

lnAr(£)c<-(-) fin -J , (30) 

(£ « T). Of course, probing experimentally the states 
with energies much smaller than the temperature is not 
an easy task, so the expression (|30[) is mainly of formal 
interest. 



The number of electron states within the tail between 
E = and E = Eq can be estimated by a diagrammatic 
approach suggested in Refl24l Section 4.3. By analogy 
with Eq.(152) of that work, one can estimate the electron 
concentration per atom as 

n c a gT/n 2 . (31) 



IV. FLUCTUATION EFFECTS IN THE CASE 
OF ELECTRONS PREFERRING CRUMPLED 
SURROUNDING 

Consider now the case of a negative coupling constant 
which physically means that corrugations result in a gain 
of electron energy. The Hamiltonian (|2|) cannot be used 
here since its spectrum is not bounded from below. In- 
stead we will use the Hamiltonian 

H = -^-V 2 + g -^(v)-^) 2 (32) 

where g is still positive and (po is the value of |V/i| 2 op- 
timal for electron hopping; the coupling constant is now 
go = — and E = is still the lowest possible en- 
ergy. Using the Hubbard-Stratonovich transformation, 
the partition function of electron in a given field tp (r) 
can be represented as (again, H = m = 1): 



p 

Ze = j Pr(r)exp -~ J 



dr(r) 
dr 



dT - g / b( r ( T )) - <Pa\ dT 



J V 1 {t) J 2?r(r)cxp -1 J 



dr(r) 

dT 



dr - i j 7 (r) [<p (r (t)) - (p ] dr - i / 7 2 (r ) ,h 



(33) 



Repeating the derivation of Eq. p^j) we will hnd the only 
difference: g m is replaced by the product i m Ylj=i 7 ( T i) 
averaged over fluctuations of the Gaussian field j(t). 

Continuing the transformations leading from Eq. (|19[) 
to Eq.(|22)) one obtains formally the same equation but 
with the replacement of g by i"fo, 



7o 



P 



7(r)dr. 



(34) 



An optimal value of this random field corresponding to 
the extremum of the exponent in Eq. (|35[) is equal to igtpo 
and the fluctuations are of the order of y/gT, that is, neg- 
ligible if T <C gifio- Thus, the effective coupling constant 
for the fluctuation contributions is just equal to .go- All 
the conclusions of the previous section about the size of 
flexuron and criterion of autolocalization remain valid. 
The mean-field expression for the shift of the band edge 
is now different. As follows directly from Eq. 



At the end, as follows from Eq. (|3"3"l) . one has to average 
the result over the Gaussian random static field 70 dis- 
tributed with the probability function 



■ exp 



2ff 



(35) 



-En 



(36) 
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V. DISCUSSION AND CONCLUSIONS 

We have shown that the electron can be self-trapped 
by bending fluctuations assuming that the parameter 
Eq. (P7)l is large enough. One can expect g to be of the 
order of the electron bandwidth W (cf. the discussions 
for graphene reviewed in RefQ). For the Young modu- 
lus a natural estimation is Y ps E co h/a 2 where E co h is 
a cohesive energy and a is a lattice constant. For co- 
valently bonded membranes (such as graphene or h-BN) 
W fa E co h- This means that the effective mass should 
be of the order of the free electron mass, and therefore 
it is hardly to expect the formation of flexuron for, e.g., 
gapped graphene with a mass m <C H 2 /Wa 2 . On the 
other hand, the criterion of autolocalization should be 
easy to fulfill for a soft membrane with E coh <§C W. 

Assuming that we are in the regime of the autolocal- 
ization, the electron injection into crystalline membranes 
at finite temperatures does not result in the conductivity, 
the states will be (auto)localized, even without external 
disorder. This will be the case till the whole fluctuation 
tail will be occupied. The number of the states in the fluc- 
tuation tail is given by Eq. (pffl) . so the transition to the 
conducting states happens at n = n c oc T. The higher 
the temperature the more electrons should be injected. 

Strictly speaking, the flexuron states are completely 
localized only if one assumes that the fluctuations are 
static. In an ideal system, the flexuron can move together 



with the anomalously corrugated region, but very slowly, 
for typical phonon times (cf. the discussion for the case 
of fluctuonii). One can expect, however, that even very 
weak extrinsic disorder will localize so slow particle. 

If the flexuron forms, its size is of the order of 1/q* 
where q* (J5J) is a crossover point from harmonic to anhar- 
monic regime, so, the bending fluctuations are stabilized 
by electrons just at the border of the strong anharmonic- 
ity. 

A very interesting issue for future studies is the case 
of finite electron (flexuron) concentrations. By anal- 
ogy with magnetic semiconductors, one can expect a 
tendency to the phase separation and electron droplet 
formation^— . Electronic mechanisms of stabilization of 
corrugations were discussed in a context of graphene^^^ 
but in a broad-gap semiconductors considered here 
physics should be essentially different. Probably, the flex- 
uron formation is a proper term in this case. 
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